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1.  Introduction 


A  study  was  recently  made  of  the  equilibrium  structure  of  non-force-free  (NFF)  mag¬ 
netic  flux  ropes  placed  in  an  ambient  medium  where  the  pressure  gradient  force  is  nonzero 
(Chen  2012,  henceforth  referred  to  as  Paper  1).  The  work  is  formulated  using  standard 
magnetohydrodynamics  (MHD)  and  develops  a  new  physical  constraint  termed  “minimum 
complexity.”  Analytic  solutions  were  given  satisfying  specified  physical  conditions  with  a 
limited  parameter  study.  The  present  report  extends  the  discussion  in  Paper  1  and  provides 
a  more  comprehensive  documentation  of  the  solution  set  as  well  as  an  additional  family 
of  solutions  that  has  not  been  previously  published.  The  basic  physics,  however,  remains 
unchanged. 

Magnetic  flux  ropes  are  important  building  blocks  of  the  solar  corona  and  interplanetary 
medium.  A  flux  rope  is  a  self-organized  magnetized  plasma  structure  consisting  of  an  electric 
current  distribution  J(x)  and  its  magnetic  held  B(x)  embedded  in  an  ambient  medium  of 
pressure  pa.  The  magnetic  held  enables  a  hux  rope  to  maintain  its  structural  identity  separate 
from  the  ambient  plasma.  In  natural  plasma  environments,  hux  ropes  reside  in  background 
plasmas.  They  are  generally  curved,  with  the  local  geometry  characterized  by  major  radius 
R  and  minor  radius  a.  For  R/a  3>  1,  the  minor  radial  forces  can  be  approximated  by  those 
acting  on  a  straight  cylinder  (e.g.,  Shafranov  1966).  Although  this  approximation  incurs  an 
error  of  0{a/R )  and  becomes  less  accurate  as  R/a  approaches  unity,  the  simple  configuration 
can  capture  the  essential  structure  of  hux  ropes. 

Within  the  framework  of  MHD  in  which  the  internal  structure  of  an  equilibrium  hux 
rope  in  the  absence  of  gravity  is  determined  by  (l/c)J  x  B  =  Vp,  where  J  =  (c/47t)V  x  B. 
A  particular  limiting  case,  the  “force-free”  (FF)  approximation  defined  by  c_1  J  x  B  =  0,  has 
been  widely  studied  because  of  mathematical  simplicity.  These  mathematical  FF  solutions 
(e.g.,  Lundquist  1950),  however,  make  no  reference  to  the  ambient  pressure  and  for  specihc 
applications,  may  not  satisfy  the  physical  boundary  conditions  dictated  by  a  magnetic  hux 
rope  embedded  in  a  plasma  medium  of  pressure  pa.  The  method  presented  in  Paper  1 
self-consistently  solves  the  basic  one-dimensional  (ID)  MHD  equation  subject  to  specihed 
asymptotic  ambient  pressure  pa.  Furthermore,  the  solutions  satisfy  requirements  that  (1) 
physical  quantities  be  everywhere  continuous  and  nonsingular,  (2)  pressure  p(r)  be  real 
and  nonnegative  (and  smoothly  converge  to  the  ambient  value  pa),  and  (3)  the  magnetic 
held  prohle  have  “minimum  complexity.”  No  specihc  radial  dependence  of  p(r)  or  B(r) 
is  prescribed.  The  derivation  of  the  method  shows  that  the  solutions  can  be  specihed  by 
two  dimensionless  parameters,  B*  =  Bt/ (&itpa)1^2  and  B*  =  Hpa/(87rpa)1/2,  where  Bt  is  the 
toroidal  (locally  axial)  held  averaged  over  the  minor  radius  and  Bpa  is  the  poloidal  (locally 
azimuthal)  held  at  the  edge  of  the  current  channel,  r  =  a.  The  requirement  of  real  and 
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nonnegative  pressure  defines  equilibrium  boundaries  beyond  which  no  physical  equilibrium 
solutions  exist.  Paper  1  also  demonstrated  that  the  average  internal  pressure  p  can  be 
separately  specified.  Some  basic  properties  of  the  self-consistent  equilibrium  were  illustrated 
with  one  solution  for  a  specified  value  of  ps  =  0.85.  One  value  of  ps,  however,  determines  only 
one  sub-family  of  solutions.  In  the  present  report,  the  treatment  is  extended  to  provide  a 
fuller  examination  of  NFF  equilibrium  solutions  in  parameter  regimes  that  were  not  discussed 
in  the  earlier  report. 

The  basic  formulation  is  given  in  Paper  1  and  is  not  repeated  here.  The  present  report  is 
not  self-contained  but  is  intended  to  supplement  Paper  1.  The  specific  equilibrium  solutions 
given  here  can  be  used  to  infer  the  qualitative  properties  of  solutions  that  are  not  specifically 
presented.  By  way  of  defining  symbols  and  physical  quantities,  an  abbreviated  derivation  of 
the  relevant  equations  is  given  in  Section  2  below.  Sections  3-5  provide  the  new  material. 


2.  One-Dimensional  Non-Force-Free  (NFF)  Flux  Ropes 

The  forces  acting  on  a  straight  cylindrical  flux  rope  are  considered  in  the  framework  of 
standard  MHD.  In  equilibrium  the  applicable  equations  are 

-J  x  B  —  Vp  =  0,  (1) 

c 

J=^VxB.  (2) 

V  ■  B  —  0.  (3) 

Gravity  is  neglected.  We  use  a  cylindrical  coordinate  system  with  the  axis  of  the  flux  rope 
aligned  with  the  z  (toroidal)  direction,  as  illustrated  in  Figure  1.  The  radial  coordinate 
measured  from  the  axis  of  the  cylinder  is  r,  with  the  electric  current  J  limited  to  r  <  a. 
The  poloiclal  (locally  azimuthal)  angle  is  denoted  by  6.  The  treatment  is  ID,  and  we  impose 
d/dz  —  d/d6  —  0.  Pressure  is  p(r)  with  the  asymptotic  ambient  pressure  pa  taken  to  be 
constant  consistent  with  the  neglect  of  gravity. 

The  current  density  has  components  Jp{r )  and  Jt(r)  in  the  poloiclal  and  toroidal  direc¬ 
tions,  producing  toroidal  Bt(r)  and  poloiclal  Bp(r )  magnetic  field  components,  respectively. 
Because  of  the  assumed  cylindrical  symmetry,  the  magnetic  field  has  the  general  structure 
B(r)  =  (0,  Bp(r ),  Bt{r ))  with  J(r)  =  (0,  Jp(r),  J*(r)),  and  the  condition  V  •  B  =  0  is  trivially 
satisfied.  Inside  the  current  channel  (r  <  a),  the  magnetic  field  “lines”  are  helical.  This 
is  schematically  illustrated  in  Figure  1  with  one  representative  helical  curve.  Outside  the 
current  channel  (r  >  a),  we  have  J(r)  =  0  and  Bp(r)  =  Bpa-{r/a)~l  and  Bt(r )  =  0  regardless 
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of  the  profile  of  J(r)  for  r  <  a.  Here,  Bpa  =  Bp{r  =  a).  By  flux  rope,  we  refer  to  the  entire 
system  consisting  of  the  current  channel  defined  by  J(r)  ^  0  (r  <  a)  and  the  magnetic 
field  including  that  for  r  >  a  attributable  to  this  J(r)  according  to  Ampere’s  law.  The 
flux  rope  is  embedded  in  a  plasma  medium  of  density  n  and  temperature  T,  with  pressure 
p(r)  =  Pi(r)  +  pe(r)  =  2 nkT,  where  we  have  assumed  Tt  =  Te.  Subscripts  i  and  e  refer  to 
ions  and  electrons,  respectively. 


In  this  cylindrical  coordinate  system,  equation  (1)  takes  on  the  form 

dp  =  _J_d&  _ 

dr  8n  dr  in  r 

where  B2  =  B2p-\-  B 2  and  d/d6  =  0  and  d/dz  =  0  have  been  used.  Now  define  the  dimen¬ 
sionless  variable  x  =  r/a  and  form  functions  bp(x)  and  bt{x)  such  that 

Bp(x )  =  Bpabp(x),  Bt{x)  =  Btbt(x).  (5) 


Here,  Bpa  =  Bp(x  =  1)  and  Bt  is  the  toroidal  field  averaged  over  the  cross-section, 

Bt  =  2n  f  Bt(r)rdr/na2  =  2Bt  f  bt{x)xdx. 

Jo  Jo 

These  definitions  lead  to  the  general  properties 

bp{  1)  =  1,  bt  =  2  [  bt(x)xdx  =  1, 


and  Ampere’s  law  takes  on  the  dimensionless  form 


dbt(x)  Id 

3-{x)  =  3t{x)  =  -xTxxb”(x)' 


(6) 

(7) 

(8) 


with  the  dimensional  current  densities  given  by  Jp(x )  =  (c/47r)(Ht/a)jp(a;),  and  Jt(x )  = 
( c/in)(Bpa/a)jt(x ).  Ampere’s  law  leads  to  Bpa  =  2It/ca,  where  It  =  2na 2  xJt(x)dx  is  the 
total  toroidal  current  integrated  over  the  cross-section  of  the  flux  rope.  This  is  equivalent  to 

i 

xjt(x)dx  =  1,  (9) 

which  is  in  turn  equivalent  to  bp(  1)  =  1.  Localization  of  the  current  to  x  <  1  with  no 
discontinuity  at  x  —  1  requires 


jp(x  >  1)  =  0,  j,(x  >  1)  =  0. 


(10) 


All  the  solutions  of  this  system  can  be  characterized  by  the  following  two  dimensionless 
parameters: 


n*  _  Bpa  *  _  Bt 

p~  (87TPa)V2’  *  -  (8npa)1/2  ' 


(11) 
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2.1.  Minimum  Complexity 


Equation  (4)  with  the  physical  conditions  still  admits  an  infinity  of  solutions.  In  order 
to  construct  unique  solutions,  the  functional  form  of  bt(x)  and  bp(x )  must  be  constrained. 
Previously,  the  FF  equation,  J  x  B  =  0,  has  been  solved  using  the  assumption  that  J(x) 
and  B(x)  have  exactly  the  same  functional  form,  satisfying  J(x)  =  aB(r)  at  every  point 
x,  where  a  is  a  constant  (Lundquist  1950).  Typically,  this  mathematical  solution  is  used 
without  reference  to  specific  boundary  conditions.  The  so-called  Lundquist  solution  may  be 
applicable  in  certain  confined  laboratory  systems  with  closed  magnetic  surfaces  by  assuming 
that  turbulent  dissipation  would  render  a  constant  in  space  while  leaving  the  total  magnetic 
hclicity  unchanged  (Taylor  1974)  or  by  assuming  that  the  system  is  in  the  lowest  magnetic 
energy  state  (Woltjer  1958).  Paper  1  adopted  an  alternative  hypothesis:  a  flux  rope  “relax” 
to  a  state  with  “minimum  complexity”  subject  to  imposed  conditions.  The  system  need  not 
be  closed,  and  no  dissipation  is  required.  Physically,  this  is  a  conjecture  motivated  by  the 
expectation  that  a  plasma  system  evolves  to  minimize  gradients  in  B(x)  and  p(x)  everywhere 
while  conforming  to  some  external  constraints  imposed  by  the  macroscopic  system. 

Mathematically,  minimum  complexity  solutions  are  constructed  as  follows.  First,  the 
cross-sectional  radius  is  characterized  by  one  length  scale,  a.  Accordingly,  bp(x )  and  bt(x ) 
are  expanded  in  terms  of  re  =  r/a  and  demand  that  the  polynomial  representation  satisfy  the 
required  constraints  with  the  smallest  number  of  terms.  This  ensures  that  there  are  fewest 
extrema.  Thus,  we  seek  solutions  of  the  form 


M 


bp{x)  =  ^  hr 


,x 


and 


m= 0 


N 


b,(x)  =  V 


(12) 


(13) 


n= 0 


and  choose  the  smallest  M  and  N  according  to  the  number  of  constraints.  From  equation  (8), 
we  obtain 


N 


jP(x)  =  -  ^  ndnxn  1, 


n= 0 


and 


M 


jt(x)  =  1)h^m  1- 


(14) 


(15) 


m= 0 


In  this  formulation,  the  minimum  complexity  hypothesis  is  a  macroscopic  constraint,  not  a 
condition  imposed  on  each  point  x  in  space. 
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3.  Alternative  Solution  1 


The  physical  constraints  listed  above  uniquely  determine  the  functional  forms  of  bt(x) 
and  bp(x)  subject  to  the  mathematical  condition  that  the  lowest-order  powers  in  x  be  used. 
There  is,  however,  a  degree  of  freedom  in  choosing  the  specific  form  of  the  functional  forms. 
In  this  section,  this  freedom  is  illustrated. 


The  family  of  solutions  SI  discussed  in  section  3.1  of  Paper  1  uses  the  lowest  three 
nonzero  terms  in  xn  for  bt(x)  and  bp(x).  Instead,  one  may  use  the  lowest  three  nonzero  even 
powers  for  bt(x)  and  nonzero  odd  powers  for  bp(x).  Accordingly,  we  choose  n  =  0,2,4  and 
m  =  1,  3,  5.  The  solution  of  the  equation  of  equation  (4)  is  then 


and 


bt  (x) 


bp(x) 


r  3  (1  —  2x2  +  a:4),  x  <  1 
k  0,  x  >  1, 

3x(l  —  x2  +  x4/3),  x  <  1, 
x_1,  x  >  1. 


(16) 


(17) 


We  will  refer  to  this  solution  as  Alternative  Solution  1,  denoted  by  Sla.  Ampere’s  law  yields 


jp(x)  =  12x(l  -x2),  (18) 

jt{x)  =  6  (1  -  2x2  +  x4),  (19) 

inside  the  current  channel  (x  <  1).  As  before,  outside  the  current  channel,  we  have  jt(x)  = 
0  and  Jp(x)  =  0  for  all  x  >  1. 

The  normalized  self-consistent  pressure  is  given  by 


p(x)  =p0-  [j Bfb2p(x )  +  B?b2t(x)  -  9 B*2]  -  9 B*2x2  (l  -  x2  +  ^x4  -  ^x6  +  -^x8 

for  inside  the  current  channel. 

Outside  the  current  channel,  the  pressure  is 

47 

p(x)  =po-  —B*p2  +  9 Bf  =  1,  x  >  1, 
which  is  manifestly  constant.  This  leads  to 


x  <  1, 
(20) 


(21) 


Po  =  1  +  (47/10)Bp*2  -  9  Bf. 


(22) 
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The  volume- averaged  internal  pressure  is  given  by  p  =  2  xp(x)dx.  For  Sla,  it  is  simply 

p  =  i  +  b;2-(Ab;2.  (23) 

Thus,  the  equilibrium  solutions  for  a  given  value  of  p  lie  on  the  curve  given  by 

/Q  V/2 

b;  =  (-B?  -  1 +p)  .  (24) 

By  demanding  B*  >  0,  this  equation  shows  that  p  >  0.8  for  all  physically  allowed  choices  of 
B*  and  B*.  Here  recall  that  for  B*  =  0,  the  allowed  solutions  have  0  <  B%  <  1/3. 

The  allowed  equilibrium  solutions  with  non-negative  pressure  satisfy  the  condition 


This  boundary  is  similar  to  curve  A  in  Figure  2  of  Paper  2  as  are  the  other  boundary  curves 
in  the  B*-Bf  space. 

Figure  2  shows  a  particular  solution  of  this  family  specified  by  B*  =  0.5  and  B*  = 
0.69189841.  This  solution  is  the  counterpart  of  Sl.l  discussed  in  Paper  1  (Figure  3,  Paper 
1)  in  that  it  lies  on  the  asymptote  to  the  equilibrium  (counterpart  of  line  B  in  Figure  2  of 
Paper  1)  with  B*  =  0.5  and  has  p(0)  =  1.  Comparing  this  solution  with  Sl.l,  it  is  clear 
that  the  field  configurations  and  current  distributions  are  virtually  identical.  This  can  also 
be  seen  from  the  value  £  =  1.7448  for  SI. la,  which  is  similar  to  £  =  1.7690  for  Sl.l,  where 

£  =  ( Bt)max/(BP )  max •  (26) 

The  most  visible  difference  is  that  for  SI. la,  the  peak  of  Bt(x )  is  slightly  broader  and  Jp(x ) 
is  peaked  slightly  farther  away  from  the  axis  than  for  Sl.l.  For  solution  SI. la,  we  obtain 
Vp  =  6.0637  x  10  2.  The  quantitative  similarity  between  Sl.l  and  SI. la  is  due  to  the  fact 
that  the  low-order  powers  x 3  and  xA  are  similar  in  the  range  0  <  x  <  1.  Figure  3  shows  the 
Alfven  velocity  Va(x)  and  plasma  /3(x)  for  SI. la,  which  resemble  those  of  Sl.l.  For  SI. la, 
we  find  0(0)  =  0.4444  and  0  =  1.2071. 

For  S2  and  S3,  similar  freedom  in  choosing  the  power  of  x  exists.  If,  however,  higher- 
order  terms  are  introduced,  additional  extrema  in  the  current  channel  may  occur,  violating 
the  minimum  complexity  condition.  Thus,  the  ability  to  choose  different  functional  forms 
for  bt(x)  and  bp(x)  is  limited. 


4.  Force-Free  Limits 


In  this  section,  we  extend  the  discussion  on  the  FF  limit  and  low-/?  limit  previously 
given  in  section  4  of  Paper  1,  where  the  limiting  properties  of  the  class  of  solutions  SI  were 
calculated  along  a  line  through  the  origin  with  a  slope  of  B*/B *  =  (350/159) 1^2  (Figure 
2,  Paper  1).  Specifically,  consider  /3(x)  and  pressure  gradient  dp(x)  j dx  averaged  over  the 
volume  of  the  current  channel: 


j3  =  2  f  x/3(x)dx,  and  Vp  =  2  f  x 


dp[x ) 


dx 


dx. 


(27) 


Here,  the  local  plasma  is  given  by  {3(x)  =  8tt p(x)  /  B2  (x) ,  and 


Bfb?{0y  and  ^  Bf  (28) 

Noting  that  (5  =  (2 u^/Va)2,  nonuniform  (3(x )  suggests  that  a  flux  rope  is  not  “monolithic” 
in  its  structural  robustness  as  it  evolves  over  time.  Here,  Vth  =  ( kT/m )1//2  is  the  thermal 
speed. 

Shown  in  Figure  4  are  the  quantitative  properties  of  /3  and  T>p  for  the  two  distinct  limits 
B*  — >  0  and  B*  oo  of  the  family  of  solutions  SI.  Figure  4a  shows  f3  and  Vp  as  well 
as  /3(0)  and  /3(1),  defined  by  equations  27  and  28,  respectively,  along  a  line  through  the 
origin  with  slope  B*/B*  =  2/3,  where  the  values  of  f3  (solid  and  dash-dot)  and  Vp  (dashed) 
are  shown  on  the  left  and  right  vertical  axes,  respectively.  On  this  line,  all  solutions  have 
£  =  3.0370  and  identical  field  profiles.  This  line  intersects  the  equilibrium  boundary  (curve 
A  in  Figure  2,  Paper  1)  at  B*  =  0.223873.  This  value  is  indicated  by  the  vertical  line  A 
in  Figure  4a,  where  the  solution  yields  p  =  0  and  therefore  /5 (0)  =  0,  while  /3  =  8.3916 
(left  axis)  and  /3(1)  =  19.9524.  This  solution  has  T>p  =  0.3638  at  A.  Beyond  this  point,  no 
physical  solutions  exist.  The  limit  B*  — >  0  leads  to  f3  — »  oo,  where  the  flux  rope  has  Vp  =  0, 
i.e.,  d,p(x)/dx  =  0  for  all  x.  This  is  the  FF  limit. 


Figure  4b  shows  the  same  quantities  along  a  line  of  slope  B*/B*  =  2.  Again,  in  the  limit 
Bp  — >  0,  we  obtain  /3  — >  oo  and  T>p  — >  0.  Flux  ropes  are  nearly  FF  and  have  high  fd  in  the 
neighborhood  of  this  limit.  According  to  equation  (28),  the  limit  B*  — >  oo  leads  to  /3(1)  =  0. 
On  the  axis  of  the  flux  rope,  however,  (3  converges  to  /5(0)  =  0.8172  with  /3  =  0.2323  (left 
axis).  The  flux  rope  becomes  highly  NFF  in  this  limit  with  T>p  — >  oo. 

In  the  solar  physics  literature,  the  common  adage  is  that  the  coronal  plasma  /3  is  so 
small  (/3  — >  0)  that  the  force-free  approximation  (c-1J  x  B  =  Vp  =  0)  is  valid.  The  above 
results  and  those  of  Paper  1  contradict  this  adage.  Because  this  adage  implicitly  underlies 
the  solar  physics  literature,  it  warrants  some  attention.  By  remote  sensing  techniques,  the 
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coronal  pressure  pa  can  be  estimated  at  various  heights,  but  the  local  coronal  magnetic  held 
cannot  be  measured  at  this  time.  Instead,  the  magnitudes  B  =  |B|  of  coronal  magnetic 
holds  are  estimated  from  photospheric  magnetograms  using  various  simplifications  such  as 
the  FF  approximation.  Using  the  estimated  “typical”  value  B'  of  the  coronal  held  and  the 
estimated  pa,  an  effective  plasma  (3  is  defined  by  f3'  =  87Tpa/B'2.  Depending  on  the  value 
of  B'  chosen  in  these  estimates,  (3'  may  indeed  be  small.  The  magnetic  held  of  a  plasma 
structure  embedded  in  a  plasma  medium,  however,  cannot  be  specihed  a  priori.  Instead, 
B  must  be  determined  self- consistently  based  on  actual  force  balance.  We  have  defined  a 
self-consistent  effective  plasma  /3  as  follows: 

at  _  8nPa 

1  ~  <  B(x )2  >’ 


where  B(x )  is  the  self- consistently  calculated  magnetic  held  in  force  balance  with  pa  and 
<  ...  >  is  the  volume  average  over  the  current  channel.  This  effective  f3'  is  plotted  in 
Figures  4a  and  4b  as  dotted  curves.  The  plots  show  that  /?'  behaves  similarly  to  the  local  (3 
and  (3  along  the  two  lines  in  the  parameter  space:  /?'—>■  oo  in  the  FF  limit  as  T>p  — y  0.  The 
usual  adage  is  no  more  valid  in  terms  of  (3 '  than  in  terms  of  the  actual  f3.  This  is  consistent 
with  the  conclusion  given  in  Section  4  of  Paper  1. 


5.  Specification  of  Average  Internal  Pressure  p 

For  S3,  the  parameter  space  structure  depends  on  the  specified  value  of  ps.  Section  3.3 
of  Paper  1  describes  in  detail  one  particular  solution  belonging  to  the  sub-family  of  solutions, 
S3.0,  defined  by  ps  =  0.85.  In  this  section,  a  more  comprehensive  discussion  of  this  family 
of  solutions  is  given  for  several  values  of  ps.  Figure  5  shows  the  BJ  B*  parameter  space 
for  S3.0.  Several  particular  solutions  are  indicated  by  the  numbered  dots,  These 

solutions  will  be  designated  S3.01-S3.05,  respectively.  The  boundary  curve  A  is  defined  by 
p(0)  =  0  with  p(0)  <  0  on  the  left  of  A.  To  the  right  of  curve  B,  pressure  p(x)  develops 
a  non- zero  imaginary  part.  The  short  line  segment  “a”  indicates  the  value  of  B*  beyond 
which  no  physical  solutions  exist,  and  the  open  diamond  represents  the  only  solution,  S3. 06, 
with  real  and  nonzero  p(x)  at  this  B*,  which  is  B*  =  1.002834.  The  line  segment  “6”  is 
at  B*  =  0.3350  below  which  all  solutions  between  A  and  B  are  physically  acceptable.  The 
open  diamond  (6)  is  at  B l  =  0.360  and  B*  =  0.3350.  Above  b,  some  solutions  near  A  and  B 
are  not  physical  because  pmin  becomes  negative. 

S3. 01  has  been  presented  in  Paper  1  and  will  not  be  discussed  here.  Figure  6  shows 
the  magnetic  field,  electric  current,  and  pressure  profiles  of  S3. 02,  where  the  constant-a 
Lundquist  solution  is  shown  (dash-dot  curves)  for  comparison.  This  solution,  indicated  by 


“2,:  in  Figure  5,  is  slightly  to  the  right  of  equilibrium  boundary  A.  For  this  solution,  pmin  ~  0 
occurs  at  x  =  0.082.  The  average  pressure  gradient  is  T>p  =  0.4075.  Figure  7  shows  the  Alfven 
speed  Va(x)  (solid  curve,  left  axis)  and  plasma  f3  (dashed  curve,  right  axis).  Because  of  the 
relatively  low  pressure  near  the  axis,  we  have  <  1.  For  this  flux  rope,  /3min  =  8.3383  x  10~3 
and  VAmax  occur  at  x  —  0.082,  and  (3  exceeds  unity  at  about  x  =  0.50,  reaching  j3  ~  5  at 
x  —  1.  Overall,  /3  =  2.5560. 

Figures  8-13 — solutions  labeled  “3”-“5”  in  Figure  5 — describe  a  variety  of  representative 
solutions  in  the  physically  allowable  parameter  space,  showing  the  profiles  of  magnetic  held, 
current,  pressure,  Ahven  speed,  and  plasma  j3.  ffere,  dash-dot  curves  show  the  magnetic 
held  components  of  the  constant-a  Lundquist  solution  for  comparison.  These  solutions  are 
chosen  to  be  near  the  boundary  curves  A  and  B. 

Above  b ,  as  boundary  A  is  approached  from  the  right,  pmm  occurs  away  from  the  axis 
so  that  pmin  <  p( 0).  This  leads  to  /3mm  becoming  negative  before  p(0)  becomes  zero,  i.e.,  at 
a  point  to  the  right  of  A.  Figures  8  provides  an  example.  This  is  in  contrast  to  the  behavior 
below  b  where  pmin  occurs  at  x  =  0  on  the  boundary  curve  A.  Above  “5,”  pmm  also  becomes 
negative  before  p(x)  becomes  complex.  Above  b ,  there  is  a  band  adjacent  to  each  of  A  and 
B  where  the  solutions  develop  pmm  <  0  and  become  unphysical.  As  B*  is  increased  between 
b  and  a,  these  bands  broaden,  and  pmin  <  0  sets  in  farther  away  from  A  or  B  until  at  b 
where  B*  =  0.3350  only  one  physical  solution  exists.  Figures  12  and  13  show  the  properties 
of  a  representative  solution  slightly  below  a.  Note  that  p(0)  increases  with  B*.  The  solution 
labeled  “a”  is  not  significantly  different  from  S3. 05  in  magnetic  field  and  pressure  profiles 
and  is  not  shown.  The  solution  at  “a”  is  designated  as  S3. 06,  for  which  B*  =  1.002834,  the 
numerical  values  of  /3mm  =  7.2004  x  10~8  and  {V a)™, ax! V aq  =  2.6435  x  103  are  obtained.  This 
shows  that  to  within  numerical  accuracy  j3min  0  and  (V7i)max  — »  oo  for  the  true  upper 
limit  of  the  physically  allowed  value  of  B*.  For  the  S3.0  sub-family  of  solutions,  there  exists 
no  FF  limit  because  ps  <  1,  and  the  limit  B*  — »  0  and  B*  — »  0  is  not  physically  allowed. 

As  the  specified  value  of  ps  is  changed,  the  parameter  space  boundaries  change,  although 
the  general  structure  remains  similar.  Shown  in  Figure  14  is  the  parameter  space  for  ps  = 
0.96395,  for  which  the  boundary  curve  A  is  a  straight  line  through  the  origin.  This  value 
has  been  numerically  determined  as  follows.  A  sequence  of  solutions  with  decreasing  Bf  and 
B*  was  evaluated,  and  for  each  pair  of  values  ps  is  adjusted  until  p(0)  <  10-6  is  obtained, 
meaning  that  the  solution  is  on  the  boundary  A.  This  sequence  was  carried  out  down  to 
B*  =  1.0  x  10-5  and  B*  =  1.0  x  10~5.  This  sub-family  of  solutions  is  referred  to  as  S3.1, 
and  specific  solutions  are  further  distinguished  by  an  additional  digit  “n”  as  in  S3.  In.  The 
boundaries  A  and  B  and  the  points  a  and  b  have  the  same  significance  as  the  corresponding 
features  in  the  parameter  space  for  the  sub-family  S3.0  (Figure  5).  For  S3.1,  solutions  above 
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B*  =  1.3601  (“a”)  are  unphysical,  while  those  below  B*  =  0.7408  (“&”)  are  all  physical. 

Figures  15-21  present  a  number  of  representative  solutions  distributed  in  the  allowed 
parameter  space  region,  corresponding  to  solutions  “l”-“4”  shown  in  Figure  14,  referred  to 
as  solutions  S3.ll-S3.14,  respectively.  For  S3. 11  (“1”  in  Figure  5),  defined  by  B*  =  0.14 
and  B*  =  0.040,  the  relatively  small  value  of  B*  yields  a  small  amplitude  Bp  component, 
small  pinch  force,  and  small  pressure  gradient,  with  T>p  =  7.4360  x  10~2.  Note  that  marked 
difference  from  the  Lundquist  solution  (dash-dot,  Figure  15)  or  SI. la.  For  S3. 11,  we  find 
£  =  7.3379,  significantly  greater  than  xil  =  1.7187  for  the  Lundquist  solution.  Because  of 
the  small  Bp  component,  the  Alfven  speed  outside  the  current  channel  is  low  and  f3  is  large 
as  shown  in  Figure  16.  For  this  solution,  (3  =  2.1165  x  102. 

Figure  17  shows  solution  S3. 12  (“2”),  defined  by  B*  =  0.068  and  B*  =  0.088,  which 
is  slightly  inside  of  boundary  A.  We  find  p(0)  =  pmin  ~  0.  Note  that  Bt{x )  is  much  more 
narrowly  peaked  than  that  of  S3. 11,  with  a  narrower  dip  in  p(x).  The  average  pressure 
gradient  is  T>p  =  0.1828.  Figure  18  shows  that  the  Alfven  speed  is  also  narrowly  peaked  at 
x  =  0.  This  solution  has  f3min  =  3.2383  x  10~2  and  f3  =  75.5251. 

Solution  S3. 13,  shown  in  Figures  19  and  20,  is  similar  to  solution  S3. 03  (Figures  8  and  9), 
and  its  physical  properties  require  no  further  discussion.  This  example  shows  that  similar 
magnetic  held  profiles  may  correspond  to  different  sub-classes  of  solutions  or  parameter 
values. 

Solution  S3. 14  (Figures  21  and  22),  indicated  by  the  open  diamond  and  “4”  in  Figure  5,  is 
the  “last”  (numerically  determined)  physically  admissible  solution.  This  solution  is  specified 
by  B*  =  0.9652  and  B*  =  1.3601.  It  has,  to  within  numerical  uncertainty,  pmin  —  0  and 
/3min  ~  0,  which  occurs  at  x  =  0.3350  where  (V^rnax/kAo  =  6.2481  x  102.  For  this  solution, 
T>p  =  1.3483  and  (3  =  0.3140,  with  (3{x)  <C  1  for  x  ^  1.  This  solution  is  similar  in  its  basic 
properties  to  S3. 05. 

Figure  23  shows  the  held  and  pressure  prohles  of  solution  S3. 15  specihed  by  B*  = 
1.0  x  10-4  and  Bf  =  1.0  x  10-4.  Panel  (a)  shows  the  held  components  (solid  curves)  and  the 
magnitude  of  the  held  |B(a;)|  (dashed).  Note  that  Bp  ~  0  so  that  |B(x)|  ~  \Bt(x)\;  they  are 
virtually  indistinguishable  except  where  Bt(x )  is  negative.  The  Lundquist  held  components 
are  shown  (dash-dot)  for  comparison.  Panel  (b)  shows  the  profile  of  the  electric  current  J(x). 
Here,  the  magnitude  J(x)  =  |J(x)|  is  shown  by  the  dashed  curve,  which  essentially  overlaps 
Jp(x )  where  it  is  positive  ( x  &  0.6).  Panel  (c)  shows  that  p  ~  0  at  x  —  0,  rapidly  increasing 
to  p  —  1  at  x  ~  0.4.  The  average  pressure  gradient  is  T>p  =  1.7958  x  10-1.  This  value 
of  T>p  <C  1  is  considerably  less  than  the  maximum  pressure  gradient  (normalized),  which  is 
the  maximum  slope  of  p{x).  This  is  because  the  bulk  of  the  dip  in  pressure  is  localized  to 


x  0.2.  In  contrast,  T>p  <C  1  is  the  volume-average  of  the  pressure  gradient  in  the  entire 
current  channel,  and  dp(x)  / dx  is  weighted  by  x  in  the  integrand. 

The  Alfven  speed  (solid  curve,  left  axis)  and  plasma  fd  (dashed  curve,  right  axis)  of 
the  system  are  shown  in  Figure  24.  In  this  structure,  fd  ~  0  near  x  —  0,  but  overall 
fd  =  4.4977  x  106  S>  1,  with  fd(  1)  =  oo.  To  within  numerical  accuracy,  p(0)  =  0  and  /3(0)  =  0. 
Because  the  temperature  is  assumed  to  be  uniform  throughout  the  system,  density  n(x)  is 
directly  proportional  to  pressure  p(x).  We  find  that  Va(x)  is  essentially  a  spike  at  x  —  0 
with  an  infinitesimal  thickness,  and  fd  — *  oo  as  the  origin  is  approached.  Thus,  even  though 
Bf  and  Bf  both  vanish,  the  structure  remains  non-force-free  with  T>p  ^  0  and  fd  — >  oo. 

The  properties  of  solutions  for  ps  >  1  have  also  been  examined.  Shown  in  Figure  25 
is  the  parameter  space  structure  for  ps  =  1.15.  As  before,  curve  A  is  the  boundary  where 
p(0)  =  0,  and  solutions  acquire  a  nonzero  imaginary  curve  to  the  right  of  curve  B.  The 
line  segment  “a”  marks  the  upper  limit  beyond  which  no  physically  admissible  solutions 
are  found,  while  below  the  segment  denoted  by  ub”  all  solutions  are  physically  acceptable, 
as  discussed  earlier.  Figures  26-33  provide  a  description  of  representative  solutions  for  this 
sub-family  of  solutions,  designated  as  “l”-“3”  and  “a,”  respectively,  in  Figure  25.  They  will 
be  referred  to  as  S3.21-S3.24,  respectively.  These  solutions  show  no  remarkable  differences 
from  those  in  the  other  sub-families  discussed  above.  As  the  parameter  values  approach  the 
boundary  A,  p( 0)  vanishes.  However,  because  ps  >  1,  the  peak  in  p(x)  is  greater,  which 
makes  the  Lorentz  force  JpBt  —  JtBp  more  negative  relative  to  the  maximum  value. 

Recall  that  for  the  families  of  solutions  SI  and  S2,  Bf  — >  0  and  Bf  — >  0  is  the  FF 
limit  with  T>p  =  0.  With  ps  ^  1,  the  pressure  gradient  cannot  be  everywhere  zero  so  that 
configurations  with  nonzero  ps  are  generally  not  FF.  For  S3.0  and  S3. 2  discussed  above, 
for  example,  the  origin  (Bf  =  0  and  Bf  =  0)  does  not  correspond  to  a  physically  allowed 
solution.  For  S3.1,  the  origin  is  on  the  boundary  A.  Thus,  p(0)  =  0.  Because  p(l)  =  1,  the 
flux  rope  necessarily  has  dp(x)/dx  ^  0  and  cannot  be  FF. 

We  now  consider  a  sub-family  where  the  origin  is  in  the  interior  of  the  physically  allowed 
region.  Choose  ps  =  0.9800.  Figure  34  shows  the  basic  parameter  space  structure.  The 
solution  S3. 31  represented  by  the  dot  at  Bf  =  0.02  and  Bf  =  0.02,  is  shown  in  Figure  35. 
Note  that  d5  in  the  definition  of  bt(x),  equation  (13),  and  therefore  bt(x)  depend  on  Bf  and 
Bf.  Unlike  solutions  in  the  families  SI  and  S2,  the  Bt  component  does  not  vanish  because 
( bt)max  increases  as  Bf  — >  0.  The  Bp(x )  and  Jt(x )  components,  however,  are  unaffected  by 
the  d5  dependence  and  vanish  with  Bf  as  before.  Here,  keep  in  mind  that  bp(x )  and  jt(x) 
are  independent  of  Bf  and  Bf.  As  a  result,  the  need  to  balance  the  forces  requires  jp(x) 
to  increase  in  value,  which  means  that  bt(x)  must  support  steeper  gradients,  leading  to  a 
narrower  dip  in  p(x)  and  a  narrower  peak  in  bt(x).  Interestingly,  as  Bf  and  Bf  are  decreased, 
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the  average  pressure  gradient  increases  only  slowly.  For  S3. 31,  we  find  T>p  =  9.6916  x  10  2. 
The  pressure  profile  develops  a  narrower  and  deeper  dip.  These  solutions  are  not  FF. 

Figure  36  shows  the  Alfven  speed  (solid,  left  axis)  and  f3  (dashed,  right  axis)  for  S3. 31. 
The  former  is  highly  peaked  near  the  axis,  and  f3(x)  is  highly  nonuniform,  being  minimum 
on  the  axis  with  the  value  /3min  —  j3(0)  —  0.8217  but  f3  3>  1  throughout  most  of  the  current 
channel,  with  (3  =  8.8781  x  102.  If  a  point  closer  to  the  origin  is  chosen,  the  field  and  current 
profiles  remain  similar  except  that  maximum  values  of  Bp  and  Jt  become  smaller  relative  to 
Bt  and  Jp,  respectively,  and  /3  becomes  much  greater. 

Note  that  the  solutions  shown  above  exhibit  similar  qualitative  properties  for  all  the 
values  of  ps  used  here.  Solutions  for  ps  =  1.20  have  also  been  examined  and  found  to  be 
similar.  We  have  been  able  to  find  physically  acceptable  solutions  down  to  ps  =  0.697793; 
at  this  value  “one”  solution  has  been  found,  having  B*  =  0.4135590  and  B*  =  0.  It  has 
Vmin  =  0  at  x  =  0.3950.  For  ps  smaller  than  this  value,  we  have  found  no  physical  solutions, 
i.e.,  non- negative  and  real  p(x).  Note  that  these  values  are  determined  numerically,  and 
there  are  inherent  uncertainties  due  to  numerical  values.  Nevertheless,  the  uncertainties  are 
insignificant  on  the  scale  of  the  plots  presented. 

For  the  sake  of  completeness,  the  parameter  space  map  for  the  family  S2  is  shown  in 
Figure  37.  Comparing  this  with  Figure  2  of  Paper  1,  it  is  clear  that  the  parameter  space 
features  are  qualitatively  similar  with  all  the  same  physical  significance. 

This  work  was  supported  by  the  Naval  Research  Laboratory  Base  Program. 


-  14 

REFERENCES 

Chen,  J.  1996,  J.  Geophys.  Res.,  101,  27499. 

Chen,  J.  2012,  Astrophys.  J.  2012,  761,  179  (Paper  1). 

Lundquist,  S.  1950,  Arkiv  Fysik  2,  361. 

Shafranov,  V.  D.  1966,  Reviews  of  Plasma  Physics,  ed.  M.  A.  Leontovich,  p.  103,  Consultants 
Bureau,  New  York. 

Taylor,  J.  B.  1974,  Phys.  Rev.  Lett.  33,  1139. 

Woltjer,  L.  1958,  Proc.  Nat.  Acad.  Sci.,  44,  489. 


This  document  was  prepared  with  the  AAS  IAIRX  macros  v5.2. 


15 


|3  »  1 


Fig.  1. —  Schematic  of  a  flux  rope  in  equilibrium  with  ambient  pressure  pa.  The  cylindrical 
coordinates  9  and  z  are  shown.  The  radial  coordinate  is  measured  from  the  axis  at  r  =  0. 
The  edge  of  the  current  channel  is  at  r  =  a.  Toroidal  (poloidal)  magnetic  field  component 
Bt  (Bp)  is  shown.  The  “flux  rope”  consists  of  the  current  channel  and  the  magnetic  field 
produced  by  the  current.  The  plasma  /3  is  typically  much  smaller  inside  than  outside  the 
flux  rope.  Reproduced  from  Paper  1. 
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SI. la  BJ  =  0.5,  Bp=  0.6919 


Fig.  2. —  Solution  SI. la.  B '*  =  0.50  and  B*  =  0.69189841  so  that  p(0)  =  1.  (a)  Magnetic 
field  components  (solid  curves)  and  the  Lundquist  solution  (dash-dot).  Thick  dashed  curve 
is  |B  (x)|.  £  =  1.7448.  (b)  Current  density  scaled  to  Jt0.  The  Lundquist  current  (dash-dot). 
Solid  dashed  curved  is  |J(x)|.  (c)  Pressure  p(x )  and  dimensionless  Lorentz  force  (dashed, 
arbitrary  units).  The  average  pressure  gradient  is  Vp  =  6.0637  x  10-2  with  pmax  =  1.0572. 


Fig.  3. —  Alfven  speed  Va(x)  (solid  curve,  left  axis)  and  f3(x)  (dashed  curve,  right  axis)  for 
solution  SI. la.  (3  —  1  is  indicated.  /3(0)  =  f3min  =  0.4444  and  (3  =  1.2071. 
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Parameter  Space 


Fig.  5. —  Parameter  space  for  Particular  Solution  3.0  with  ps  =  0.85.  Solution  S3. 01  shown 
in  Figure  12  of  Paper  1  is  marked  by  “1.”  Line  segment  “a”  is  at  B*  =  1.002834,  above 
which  no  physical  solutions  can  be  found,  and  line  segment  “6”  is  at  B*  =  0.3350,  below 
which  all  solutions  are  physical.  Solutions  S3. 02-3. 05  are  indicated.  The  open  diamond  at 
“a”  is  solution  S3. 06  defined  by  B =  0.761601  B*  =  1.002834. 
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S3. 02  B(  =  0.40,  Bp=  0.445 


0.0  0.5  1.0  1.5  2.0 


X 

Fig.  6. —  Solution  S3. 02  (“2”  in  Figure  5)  given  by  B*  =  0.40  and  B*  =  0.4455.  As  required, 
p  =  ps  =  0.85.  Pmin  =  1.6439  x  10-2  at  x  =  0.082.  Vp  =  0.4075.  (a)  Magnetic  field  (solid) 
scaled  to  Bt0,  compared  the  Lundquist  solution  (dash-dot),  (b)  Current  density  scaled  to 
Jt0.  (c)  Pressure  and  normalized  Lorentz  force  (dashed).  See  Paper  1  for  S3. 01. 


Fig.  7. —  Alfven  speed  Va(x )  (solid  curve,  left  axis)  and  f3(x)  (dashed  curve,  right  axis)  for 
solution  S3. 02.  (3  =  1  is  indicated  (thin  dashed  line).  Va(x)  is  peaked  at  x  =  0.080  where 
(3  =  /3min  =  8.3383  x  10-3  and  and  ft  =  2.5560. 
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S3. 03  Bj  =  0.6,  Bp=  0.7800 


Fig.  8. —  Solution  S3. 03  (“3”)  given  by  B*  =  0.60  and  B*  =  0.7800.  Jp(x)  is  slightly 
negative  near  x  =  0  so  that  Bt(x)  is  peaked  off  axis.  pmin  =  1.0529  x  10~2  at  x  —  0.2200 
and  pmax  =  1.0962  at  x  =  0.6900.  Vp  =  5.9044  x  lO^1. 


Fig.  9. —  Alfven  speed  Va(x)  (solid,  left  axis)  and  {3(x)  (dashed,  right  axis)  for  solution 
S3. 03.  Va(x)  is  sharply  peaked  at  x  =  0.2200.  /3mm  =  2.8046  x  10~3  but  ft  =  8.6924  x  10-1. 
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S3. 04  B(  =  0.64,  Bp=  0.780 
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Fig.  10. —  Solution  S3. 04  (“4”)  given  by  B*  =  0.64  and  B*  =  0.7800.  Jp(x )  is  strongly 
negative  near  the  axis  with  a  significant  peak  in  Bt(x).  Pressure  p(x)  is  peaked  on  the  axis 
with  pmax  =  2.9268  and  /3mm  =  p(0.4)  =  1.1037  x  10_1.  Vp  =  1.0801. 


Fig.  11. —  Alfven  speed  Va(x)  (solid,  left  axis)  and  /3(x)  (dashed,  right  axis)  for  solution 
S3. 04.  Va(x)  is  peaked  at  x  —  0.40  where  /3(x)  has  a  minimum,  f3min  =  3.5829  x  10-2. 
/3(0)  =  2.5581  and  ]}  =  8.3177  x  10”1. 


Fig.  12. —  Solution  S3. 05  (“5”)  given  by  B*  =  0.73  and  B*  =  0.950.  Similar  in  character  to 
S3. 04.  This  solution  has  pmax  =  2.1130  and  Vp  =  0.9490. 
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Fig.  13. —  Alfven  speed  Va(x)  (solid,  left  axis)  and  {3(x)  (dashed,  right  axis)  for  solution 
S3. 05.  ]}  =  0.5765  with  pmin  =  1.0794  x  10“2  at  x  =  0.3400  and  /3(1)  =  1.1080. 
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Parameter  Space 


Fig.  14. —  Parameter  space  for  family  S3.1  with  ps  =  0.96395.  The  equilibrium  boundary  A 
is  straight  and  goes  through  the  origin.  Diamond  “a”  is  at  B*t  =  0.9652  and  B*  =  1.3601, 
and  line  “6”  is  at  B*  =  0.7408.  Above  a,  no  physical  solutions  exist.  Below  b  all  solutions 
between  A  and  B  are  physical. 
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S3. 11  B{  =  0.14,  Bp=  0.040 


Fig.  15. —  Solution  S3. 11  (“1”  in  Figure  14)  with  B*  =  0.14  and  B*  =  0.040.  Because 
B*J B*t  <C  1,  the  Bt  component  dominates,  with  a  peak  off  axis  and  negative  Jp  near  the 
axis.  Vp  =  7.4360  x  10~2.  The  Lundquist  solution  (dash-dot)  is  shown  in  panel  (a). 


Fig.  16. —  Alfven  speed  Va(x )  (solid,  left  axis)  and  (3(x )  (dashed,  right  axis)  for  solution 
S3. 11.  ^  >  1  throughout  the  structure  with  /3(0)  =  9.5710,  /3min  =  6.1916,  and  f3  = 
2.1165  x  102. 
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S3. 12  B',=  0.068,  Bp=  0.088 


Fig.  17. —  Solution  S3. 12  (‘2”)  with  B*  =  0.068  and  B*  =  0.088.  Both  Bt(x )  and  Jp(x )  have 
negative  segments  near  x  —  1.  T>p  =  0.1828.  The  Lundquist  field  is  shown  in  panel  (a)  for 
comparison. 


Fig.  18. —  Alfven  speed  Va(x)  (left  axis)  and  f3(x)  (right  axis)  for  solution  S3. 12.  The  Alfven 
speed  is  sharply  peaked  on  the  axis,  and  (3  rapidly  increases  from  [3mtn  =  3.2383  x  10”2  at 
x  =  0  and  becomes  much  greater  than  unity  at  x  1.  /3  —  75.5251. 
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Fig.  19. —  Solution  S3. 13  (“3”)  defined  by  B*  =  0.70  and  B*  =  1.013.  pmax  =  1.2666  at 
x  =  0.6200,  Pmin  =  5.2032  x  10~2  at  x  =  0.1800,  and  Vp  =  0.7023.  The  field  profiles  are 
similar  to  those  of  the  Lundquist  solution  (dash-dot  curves). 
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Fig.  20. —  Alfven  speed  Va(x )  (left  axis)  and  /3(x)  (right  axis)  for  solution  S3. 13. 
(V a) max /v AQ  =  2.2569.  f3  <  1  for  x  &  0.3  with  /3(0)  =  4.5160  x  10"2  and  ]}  =  5.7220  x  10"1. 


27 


S3. 14  B,'=  0.9652,  Bp=  1.3601 


Fig.  21. —  Solution  S3. 14  (“4”)  given  by  B%  =  0.9652  and  B*  =  1.3601,  denoted  by  the  open 
diamond  “a”  in  Figure  14.  The  only  physical  solution  at  B*  =  1.3601.  Jp(x)  is  negative 
slightly  off  the  axis.  T>v  =  1.3483. 


Fig.  22. —  Alfven  speed  Va(x)  (solid  curve,  left  axis)  and  {3(x)  (dashed  curve,  right  axis) 
for  solution  S3. 14.  (Vd)max/Kio  =  1-4134  x  102.  The  flux  rope  has  (3  =  3.1395  x  10-1  with 
Pmin  =  2.0231  x  10-5  at  x  —  0.34.  For  this  flux  rope,  /l  <  1  for  all  x  <  1. 


S3. 15  B*t  =1.0x1 0  4,  Bp  =  1.0x1 0  4 


Fig.  23. —  Solution  S3. 15  given  by  B*  =  10  4  and  B*  =  10  4.  The  average  pressure  gradient 
is  Vp  =  1.7958  x  10”1.  To  within  the  numerical  accuracy,  Va( 0)  —  oo  because  p(0)  ~  0. 


Fig.  24. —  Alfven  speed  Va(x)  (solid  curve,  left  axis)  and  {3(x)  (dashed  curve,  right  axis) 
for  solution  S3. 15.  Va(x)  is  essentially  a  spike  at  x  —  0,  decreasing  from  Va/Vao  =  1  to 
Va/Vao  =  0.1  at  x  ~  0.04  and  to  Va/Vao  =  0.01  at  x  ~  0.25.  /3(0)  =  6.6474  x  10-4  but 
overall  f3  =  4.4977  x  106  S>  1. 
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Parameter  Space 


Fig.  25. —  Parameter  space  for  family  S3. 2  with  ps  =  1.1500.  Line  segment  “a”  is  at 
B*  =  1.80004  and  “5”  is  at  B*  =  1.11306.  At  a,  one  physical  solution  (open  diamond) 
exist,  and  above  “a”  no  physical  solutions  exist.  Below  b  all  solutions  between  A  and  B  are 
physical.  To  the  left  of  A,  p(0)  <  0,  and  to  the  right  of  B ,  p(x)  has  a  nonzero  imaginary 
part. 
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S3. 21  Bf  =  0.12,  Bp=  0.500 
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Fig.  26. —  Solution  S3. 21  (“1”  in  Figure  25)  given  by  B*  =  0.12  and  B*  =  0.50.  Jp(x)  has  a 
negative  segment  near  x  —  1.  pmm  =  0.1466  and  pmax  =  1-5527  at  x  —  0.36.  Vp  =  0.5331. 
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Fig.  27. —  Alfven  speed  Va(x)  (left  axis)  and  (3(x )  (right  axis)  for  solution  S3. 21.  Va(x)  is 
sharply  peaked  at  x  —  0.  Except  near  the  axis,  f3  S>  1  with  f3  =  3.1029. 
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S3. 22  Bf  =  0.30,  Bp=  0.586 


Fig.  28.—  Solution  S3.22  (“2”)  given  by  B*  =  0.30  and  B*p  =  0.586.  Vp  =  0.2771. 


Fig.  29. —  Alfven  speed  Va(x)  (left  axis)  and  /3(x)  (right  axis)  for  solution  S3. 22.  Va(x)  has 
a  broad  and  gentle  peaked.  (3  >  1  throughout  the  flux  rope,  having  /3  =  1.9670. 
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Fig.  30. —  Solution  S3. 23  (“3”)  given  by  B*  =  0.80  and  B*  =  1.230.  Vp  =  0.7509. 
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Fig.  31. —  Alfven  speed  Va(x)  (left  axis)  and  /3(x)  (right  axis)  for  solution  S3. 23.  f3  <  1 
throughout  the  current  channel  with  /3  =  0.4325. 
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Fig.  32. —  Solution  shown  by  an  open  diamond  at  “a”  at  B*  =  1.22824  and  B*  =  1.79972, 
the  only  physical  solution  at  this  value  of  B*.  pmin  —  0  at  x  —  0.33.  Vp  =  1.9912. 
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Fig.  33. —  Va(x)  (solid,  left  axis)  and  {3(x)  (dashed,  right  axis)  for  solution  “a”. 
{V a) max /V ao  =  1.6869  x  103,  Pmin  —  0  at  x  —  0.33,  P  =  0.2060,  and  P{x)  <C  1  for  all 

x  <  1. 
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Parameter  Space 


Fig.  34. —  Parameter  space  for  family  S3. 3  with  ps  =  0.9800.  The  origin,  B* 
B*  =  0,  is  in  the  physically  allowed  region. 


=  0  and 
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S3. 31  Bt*  =  0.02,  Bp=  0.020 
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Fig.  35. —  Solution  S3. 31  given  by  B*t  =  0.02  and  B*  =  0.02  and  indicated  by  the  dot  in 
Figure  34.  Vp  =  9.6916  x  10“2. 


Fig.  36. —  Va(x)  (left  axis,  solid  curve)  and  {3(x)  (right  axis,  dashed  curve)  for  solution  “a”. 
Pmin  =  /3(0)  =  0.8217  at  x  =  0,  but  ]}  «  8.8781  x  102  >  1. 


Fig.  37. —  Parameter  space  for  Particular  Solution  2  to  be  compared  with  that  for  Solution 
1  shown  in  Paper  1,  where  solution  S2.1  indicated  by  point  “1”  is  discussed. 


